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If dark energy interacts with dark matter, there is a change in the background evolution of the 
universe, since the dark matter density no longer evolves as a -3 . In addition, the non-gravitational 
interaction affects the growth of structure. In principle, these changes allow us to detect and 
constrain an interaction in the dark sector. Here we investigate the growth factor and the weak 
lensing signal for a new class of interacting dark energy models. In these models, the interaction 
generalises the simple cases where one dark fluid decays into the other. In order to calculate the 
effect on structure formation, we perform a careful analysis of the perturbed interaction and its 
effect on peculiar velocities. Assuming a normalization to today's values of dark matter density and 
overdensity, the signal of the interaction is an enhancement (suppression) of both the growth factor 
and the lensing power, when the energy transfer in the background is from dark matter to dark 
energy (dark energy to dark matter). 



I. INTRODUCTION 

Advances in cosmological observations have led to 
strong evidence for non-baryonic cold dark matter and for 
a late-time acceleration of the universe, possibly driven 
by a dark energy field (see, e.g. [H, H, Q)- Dark energy 
and dark matter are the dominant sources in the 'stan- 
dard' model for the evolution of the universe. Both are 
considered essential missing pieces in the cosmic puzzle 
- and both are currently only detected via their gravita- 
tional effects. There could therefore be an interaction be- 
tween them without violating current observational con- 
straints. Furthermore, such an interaction could alleviate 
the 'coincidence' problem (why are the energy densities 
in the two components of the same order of magnitude 
today?). And interacting dark energy, by exerting a non- 
gravitational 'drag' on dark matter, introduces new fea- 
tures to structure formation, including possibly a new 
bias [1] and a violation by dark matter of the weak equiv- 
alence principle [1, Q . 

In order to pursue the possibility of interacting dark 
energy, we need to compute the effect of the interaction 
on the background expansion and on structure forma- 
tion, and then to confront the results with data. The 
main problem is that there is no fundamental theory to 
guide us as to the form of an interaction. This problem is 
in fact subsidiary to a bigger problem - that there is cur- 
rently no fundamental theory to underpin any model of 
dark energy, including non-interacting dark energy. The 
situation is somewhat similar to that in reheating of the 
universe after inflation: there is no fundamental theory 
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for the inflaton field, and no fundamental theory to guide 
us as to the form of the interactions between the inflaton 
and other fields during reheating. In this situation, wc 
are forced to adopt phenomenological models to explore 
and narrow down the space of possibilities. 

Various interacting dark ener gy m odels have been put 
forward (see, e.g. M, j j, [ j , [ToL llll . [ill H3 , Hjl EE EE 

Ezl, EI El HE Hi MM HI Ha Hi HIW)- ah f 

these models are phenomenological. Some of them are 
constructed specifically for mathematical simplicity - for 
example, models in which the energy exchange rate is 
proportional to the Hubble rate. We consider models 
which arc similar to simple models of reheating and of 
curvaton decay, i.e., where the energy exchange is in the 
form of a decay of one species into another: the energy 
exchange is linear in the energy density of the decaying 
species, and the decay rate is constant. Such models 
were introduced in [tH . |22| . where dark matter decays 
to dark energy. Here we follow [2!| and generalize the 
energy exchange to a linear combination of dark sector 
energy densities, including as special cases the decay of 
dark matter to dark energy, and the decay of dark energy 
to dark matter. 

We analyze the growth rate of structure and the weak 
lensing convergence, both of which are sensitive to in- 
teracting dark energy. In a previous work [30| . we con- 
sidered the case of dark matter decaying to dark energy. 
Here we generalize to include the case of dark energy 
decaying to dark matter. In [3(| we did not analyze 
the peculiar velocity of dark matter, but rather deduced 
the density perturbation evolution using qualitative ar- 
guments. Here we perform a a careful analysis of the pe- 
culiar velocity, which confirms the qualitative arguments 
of [3(j - but also allows us to deal with the new case of 
dark energy decay, where the peculiar velocity analysis 
is essential for deriving the correct density perturbation 
evolution. 
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Following the typical approach in non-interacting dark 
energy models, we use the parametrization of the equa- 
tion of state for the dark energy [^3, H2] , 



w(a) = w + w a (l — a) 



(1) 



In the background, a general coupling can be described 
by the continuity equations of cold dark matter (c) and 
dark energy (x), 



Pc - 
■3H(1 



3Hp c 
-w)p x 



Q, 



-Q c 



(2) 
(3) 



where w = P x /p x , the background metric is ds 2 = —dt 2 + 
a 2 dx 2 , and Q c (Q x ) is the rate of energy density transfer 
to dark matter (dark energy). Therefore Q c < (> 
0) implies that the direction of energy transfer is dark 
matter — > dark energy (dark energy — > dark matter). 

In order to avoid stringent "fifth force" constraints, we 
assume that baryons (b) and photons (7) are not coupled 
to dark energy: Qb = = Q 1 . We also neglect radiation 
since we are focusing on structure formation in the late 
universe. The energy balance equation for fluid A is 

p A + 3i?(l + w a )pa = Qa , (4) 

with w c = Wb = 0, w x = w and Q c — —Q x 5^ = Qb- 
The Friedmann equation is 



2 SttG . _ _ _ . 
H = ~^-{pc + Px + Pb)- 



Background dynamics 



(5) 



Once a form of Q c is given, the background dynamics 
are fully determined by Eqs. Q and (0, with w given 
by Eq. (JTJ) . Here we use the new model [29| 



Qc = ~ (r c p c + T x p x ) 



(0) 



where are constant energy density transfer rates. The 
special cases of Eq. ^ are: decay of dark matter to dark 
energy, i.e., T c > — T x , and decay of dark energy to 
dark matter, i.e., < = T c . Models of the pure-decay 
type have been used in reheating after inflation [33j , and 
to describe the decay of dark matter into radiation O 
or of a curvaton field into radiation 
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A complete dynamical systems analysis for Eq. ([6]) is 
given in [29| for the case w = const, i.e., w a = in Eq. (JTJ) . 
In this paper we will consider variable w, and the partic- 
ular cases when one of the transfer rates T a is zero. The 
special case previously considered in [22] corresponds to 
T x = 0, with w a — 0. As shown in [22|, when T c > 
and w a — 0, the dark energy density becomes negative 
at early times. The source of this problem is the rigidity 
of the assumption that w is constant - the problem does 
not arise for quintessence models. It can also be avoided 
with variable w for suitable choices of wq, w a in Eq. (fT]), 
as shown in 1361] - 



The case T c < (with w a = 0) avoids negative dark 
energy density, but there is no attractor solution; by con- 
trast, there is a late-time attractor when T c > [2^ ]. 
Equation @ has an exact solution for = 0, 



p c = Pcoa exp[-T c (t-t )], 



(7) 



which shows that the dark matter density is always pos- 
itive, regardless of the sign of r c . For the special case 
r c = = w a , the dark energy density is always positive, 
as can be seen from the exact solution of Eq. ([3]): 



Pxoa 



-3(l+tu) 



enp[T x (t - t 



(8) 



II. DENSITY AND VELOCITY 
PERTURBATION EQUATIONS 

Since we are interested in the late universe we can ne- 
glect anisotropic stress, and scalar perturbations of the 
flat metric are given, in Newtonian gauge, by 

ds 2 = -(1 + 2(f))dt 2 +a 2 (l -2<j))dx 2 . (9) 

The A-fluid four-velocity is 

u t X=(l-cj>,d i v A ). (10) 

Choosing the energy frame for the total four- velocity u M , 
the total velocity potential v is defined by 



PA 



PA 



Pa) 



VA 



(11) 



The covariant form of energy-momentum transfer is 



V„T- 



v A 



(12) 



The energy-momentum transfer four-vector can be split 
relative to the total four-velocity as [22| 

Q\ = Qau» +F*X, Q a = Qa + SQa, = 0, (13) 

where is the momentum density transfer rate, relative 
to it". Then it follows that F% = (0,d l f A ), where f A is 
a momentum transfer potential, and 

Qo = - [Qa (1 + 0) + SQ A ] , (14) 
Qf = a 2 di(f A + Q A v). (15) 

The general evolution equations for the dimensionless 
density perturbations S A — Spa/ Pa and the velocity per- 
turbations va are [22| : 

k 2 

5 A + 3Hc 2 sA 5 A - (1 + w A ) — v A 

a 

- 3H[3H(1 + w A ){c 2 . A ~ w a) + w a ]va ~ 3(1 + w A )(t) 
SQa , Qa r± r n _ TTI 2 



PA PA 

v A + H(l - 3c 2 A )v A 
1 



S A - 3aH(c sA - w A )v A ] 

„2 ± 



(16) 



SA 



, , < S A + ~ 

a(l + wa) a 



(1 + Wa)Pi 



{Q A [ v -(l + c 2 sA )v A ]+f A }. (17) 
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The relativistic Poisson equation is: 

*i = -3H j> - ZH 2 4> ~ 4nG (p c S c + p b S b + p x 5 x ) . (18) 

Structure formation takes place in the Newtonian 
regime, on spatial scales much smaller than the horizon 
radius, a/k <C H^ 1 . Then the gravitational potential 
and its time derivative may be neglected relative to mat- 
ter density fluctuations. This allows us to discard the <j> 
terms in Eq. (p~6|) and the first two terms on the right 
of Eq. (JTHJ). Dark energy fluctuations may be neglected 
since dark energy has a high sound speed, and does not 
cluster on sub-Hubble scales (we take the sound speed of 
dark energy to be that of a standard scalar field model, 
c 2 x = 1). In the Newtonian regime, the evolution equa- 
tions for dark matter become: 

So-— Vc = -^(5Qc~QcSc) , (19) 
a p c 

v c + Hv c + -cf> = — [Q c (v - v c ) + fc] , (20) 
a p c 

while for baryons: 

k 2 

S b v b = 0, (21) 

a 

v b + Hv h + -cf) = 0. (22) 
a 

The Poisson equation in the Newtonian limit and neglect- 
ing dark energy clustering is 

k 2 

-*<l> = -4irG (p c S c + p b S b ) . (23) 

Dark sector energy-momentum transfer 

In order to analyze structure formation in interacting 
dark energy models, we need to specify a covariant form 
of the transfer four-vectors, Q£ — —Q x , which recovers 
the energy transfer model of Eq. (J6j) in the background. 

Firstly, following and generalizing [22j , we can promote 
Eq. ([6]) to the perturbed universe directly via 

Q c := Q c + 6Q C = -T c p e (l + 5 e ) - T x p x (l + 5 X ). (24) 

However, the background model provides no guide for 
specifying the momentum transfer (which vanishes in the 
background) . It is crucial to provide a covariant prescrip- 
tion that fixes the momentum transfer, a point that is not 
always clearly recognised in the literature. Following [22| . 
we specify the momentum transfer by the covariant phys- 
ical requirement that there is no momentum transfer in 
the dark matter frame, i.e., 

Qc 1 = Qc< • (25) 

Then Eqs. (f2"4")l and (f2"S"]) completely describe the interac- 
tion. 



Regardless of the specific form of Q c , it follows from 
Eq. ([55]). using Eqs. dTJ) and JT5J, that 

fc = Qc {Vc - v) , (26) 
so that Eq. (|2"0|) becomes 

v c + Hv c + - = . (27) 
a 

This means that the dark matter velocity is the same 
as the baryon velocity and is not directly affected by 
the interaction with dark energy. In particular, the form 
of energy-momentum transfer in Eq. |25j) , for any Q c , 
ensures that there is no violation of the weak equivalence 
principle for dark matter, as shown in This will not 
be true for interaction models where Q£ is not parallel 
to <. 

III. GROWTH FACTOR AND WEAK LENSING: 
THE CASE T,; = 

For this model, 

= -r c p c (l + 5 c )u» . (28) 

It follows from Eqs. ^ and © that in the background, 
T c > corresponds to the decay of dark matter into dark 
energy, while r c < describes energy transfer from dark 
energy to dark matter. 

The equations for density and velocity perturbations 
and the Poisson equation are 

5c~— - 0, (29) 
a 

v c + Hv c + ^ = 0, (30) 
a 

k 2 (j) = -AKGa 2 {p c 8c + p b 5 b ). (31) 

From these equations, the second-order evolution equa- 
tion for 5 C is 

6 C + 2H5 C - 4nG(p c 6 c + p b 6 b ) = . (32) 

These equations have precisely the same form as those 
for the non-interacting case. However, the solutions 5 C 
are different because the background terms H and p c 
evolve differently. The interaction will produce a signal 
in the growth of structure, as measured by the growth 
factor and weak lensing measurements. Equation (|32p 
confirms the qualitative arguments used in [3fJ, where 
the velocity perturbations were not analyzed. In this 
section, we extend the analysis of [30] to include the case 

r c <o. 

For our computations, the initial conditions are set at 
the present epoch o,q = 1, and we integrate backwards 
to a = 10~ 2 . We normalize the dark matter background 
densities and density perturbations in the interacting and 
non-interacting cases to today's values, f£ c o and 5 c o- We 
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take r c = ±0.3iJo, and use an equation of state (Q} with 
w = —0.99 and w a = 0.8. These values are known [36] 
to avoid p x < in the past when r c > 0. For simplicity, 
and since we are only illustrating the effects of interact- 
ing dark energy, rather than making accurate parameter 
estimations, we will neglect the baryons in our computa- 
tions. 

In Fig. Q] we plot the energy densities Q,a = 
%-kGpaI^U' 1 . Models with dark matter decay (T c > 0) 
have higher dark matter density in the past relative to 
the non-interacting models (recall that we normalize the 
densities in the two models to the same values at cto = 1). 
By contrast, models with energy transfer from dark en- 
ergy to dark matter (r c < 0) exhibit lower dark matter 
density in the past. 

Figure [5] shows the growth function D + = 5 c /5 c q in 
the interacting models compared to the non-interacting 
case, normalized to today's value. The quantity D+ j a 
reflects the time evolution of the gravitational poten- 
tial <f>. The model with decaying dark matter, T c > 0, 
shows an enhancement of structure growth relative to the 
non-interacting model, since a lower initial value at early 
times leads to the same final value today, S c q, as the non- 
interacting model. By contrast, models with r c < show 
a suppression of structure growth. These features are 
consistent with the background evolution shown in Fig.[T] 
- recall that the difference in the growth function is de- 
termined entirely by the different background evolution, 
since Eq. (|3"2)) has the same form as the non-interacting 
version. When r c > 0, there is more dark matter in the 
past, and this leads to an enhancement in the growth 
of structure. The reverse holds for T c < 0. Note that 
the enhancement / suppression is specific to the circum- 
stance that the interacting and non-interacting models 
are normalized to have the same parameters today, i.e., 
flco and S c0 . 

The weak lensing convergence spectrum is given by (30j 



C k (£) 



9 

4c? 



d X G 2 (aH) 4 n 2 c D 2 + P(k = £/ x ), (33) 



where x( a ) — Ja da/a 2 H(a) is the comoving distance, 
P(k) is the dark matter power spectrum, and 



G( X ) 



ax PK z )-n T~ 

dx X 



(34) 



withp(z) giving the redshift distribution of lensing galax- 
ies. 

Figure [3] illustrates the impact of interacting dark en- 
ergy on the weak lensing convergence power spectra. 
The lensing power in models with decaying dark matter 
(r c > 0) is higher since the dark matter density in the 
past was higher (see Fig.[T]), compared to non- interacting 
models - assuming that the models are normalized to 
have the same values today of £l c o- The higher mat- 
ter density leads to stronger gravitational potential and 
hence a stronger light deflection. When r c < 0, the re- 
verse situation applies. 



Our results in the case r c > are consistent with those 
of [30j , where in addition the weak lensing bispectrum is 
computed. The case r c < was not considered in [30l ]. 



IV. GROWTH FACTOR AND WEAK LENSING: 
THE CASE r c = 

Weak lensing has not previously been analyzed in this 
model, for which 



-T x p x (l + 5 x )u c l . 



(35) 



It follows from Eqs. (J2J) and © that in the background, 
< corresponds to the decay of dark energy into dark 
matter, while > describes energy transfer from dark 
matter to dark energy. This model does not present the 
problem of negative dark energy density in the past. 

The velocity perturbation equation (|27[) is the same 
as in the previous model, but the density perturbation 
equation l|19p has a source term: 



Sr.- 



T P x X 

Pc 



(36) 



which will generate a linear bias between dark matter and 
baryons. This source term leads to a modified evolution 
equation for S c : 



2//(l-^gH,.-4,T^ / .,. 



1 + 

da H p c 



3w 



Tx 
H 



Px_ 

Pc 



PbSb 



0. (37) 



The evolution equation for dark matter density pertur- 
bations is modified relative to the non-interacting models 
in 3 ways: 

(1) due to the different background evolution (affecting 
the terms p c ,p x ,H); 

(2) due to the modified Hubble friction term, H — » 
H(l - T x p x /Hp c ); 

(3) due to the modified source term, which gives rise to 
a modified effective Newton constant G e ff for dark 
matter: 



2 - 3to + 1 



F, 
H 



G 3a H p c 

so that 

S c +2h(i - ^yjS c -A7TG eS p c 8 c -AnGp b 8 b = 0. (39) 

Compared to the previous model, where only modifica- 
tion (1) operates, this model shows a more complicated 
change from the non-interacting case. 
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FIG. 1: The densities of dark matter, Sl c (thin lines), and dark energy, Q x (thick lines), in the interacting models (dashed-dotted 
lines), with T c = ±0.3-ffo, and non-interacting models (solid lines), normalized to today's values. 
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FIG. 2: Linear growth function D+ = 8 c /S c o, normalized to today's value, relative to its value in a pure-matter model (D+ = a). 
The interacting models (dashed-dotted lines), with F c = ±0.3i/o, are shown in comparison to non-interacting models (solid 
lines). 



Figure [4] shows the energy densities VLa- (We again 
neglect baryons in the plots.) The model with energy 
transfer from dark matter to dark energy (T x > 0), has 
higher dark matter density in the past. The model with 
dark energy decay to dark matter (r^ < 0) has lower dark 
matter density in the past. (We normalize to today's 
value of r2 c rj.) 

In Fig. [5] we show the growth function obtained by 
solving Eq. ([55]) (neglecting baryons). As in the previous 
model, we see that an energy transfer from dark matter 
to dark energy, r x > 0, leads to enhancement of growth, 
since less initial power is needed to obtain the same final 
value as the non-interacting model. The reverse holds for 



the model with T x < 0, which suppresses growth for the 
same final value <5 c rj. These background effects are rein- 
forced by the modifications (2) and (3) to the 5 C evolution 
equation Eq. ([39f : 



• When T x > 0, the Hubble friction term is reduced, 
and the source term is enhanced, both contributing 
to an increase in 8 C relative to the non-interacting 
case. 

• For T x < 0, the friction term is enhanced and the 
source term is reduced, both contributing to a re- 
duction in 6 C . 
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FIG. 3: Weak lensing convergence power spectra £(£ + l)Ck{l) /2n in the interacting models (dashed-dotted lines), with 
T c = ±0.3-Ho, in comparison to non-interacting models (solid lines). 
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scale factor a scale factor a 

FIG. 4: The densities of dark matter, Sl c (thin lines), and dark energy, Q x (thick lines), in the interacting models (dashed-dotted 
lines), with F x — +0.3Hq, —0.2Hq, and non- interacting models (solid lines), normalized to today's values. 



Using the growth function in Eq. (|33|). we compute 
the impact of interacting dark energy on the weak lens- 
ing convergence power spectra. The results are shown in 
Fig. [6j given the normalization to today's values of f2 c o 
and S C Q, the model with > (T x < 0) exhibits en- 
hanced (suppressed) power relative to the non-interacting 
model. The cause of this is both the background ef- 
fect, i.e., since the matter density in the past was higher 
(lower) than the non-interacting case, and the additional 
reinforcing effects of modified Hubble friction and mod- 
ified Newton constant. These effects work together to 
produce stronger (weaker) gravitational potentials and 
hence a stronger (weaker) light deflection. 



V. CONCLUSIONS 

For the interacting dark energy models described by 
the covariant interaction four-vector 

QH = Qc<, (40) 

Qc = -T c p c (l + 5 C )-T x p x (l + 5 X ), (41) 

we investigated the effect of the interaction on the back- 
ground evolution of densities, on the growth factor in 
structure formation, and on the weak lensing convergence 
power spectrum. 

We considered the two cases, (I) F^ = and (II) T c = 
separately. In the general case where both transfer 
rates are nonzero, the results are a linear superposition 
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FIG. 5: Linear growth function D+ = S c /S c o, normalized to today's value, relative to its value in a pure-matter model (D+ = a). 
The interacting models (dashed-dotted lines), with F x — +0.3_Ho, — 0.2Ho, are shown in comparison to non-interacting models 
(solid lines). 
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FIG. 6: Weak lensing convergence power spectra £(£ + l)Ck(£) /2n in the interacting models (dashed-dotted lines), with 
T x — +0.3Ho, ~0.2Ho, in comparison to non-interacting models (solid lines). 



of the two cases, and will depend on the signs and relative 
magnitudes of Ta [the density perturbations will be given 
by Eq. O]. 

Case I does not produce changes in the evolution equa- 
tion for the dark matter density perturbations, so that 
the change in the growth factor and weak lensing is due 
purely to the different background evolution. In case II, 
the 5 C evolution equation has two additional modifica- 
tions: a modified Hubble friction term, and a modified 
source term, expressed as a modified effective Newton 
constant for dark matter: 

S c +2h(i - ^r^j6 c -^7TG cS p c 5 c -4TTGp b S b = 0, (42) 



2 Y 



3w 



where 

G 3a H p c 

These modifications act in the same direction as the mod- 
ification to the background evolution. 

Both cases present the same qualitative feature in the 
background evolution and in the growth of structure, as- 
suming a normalization to today's values of f2 c o and S c o'. 
for energy transfer from dark matter to dark energy, i.e. 

> 0, the dark matter density is higher in the past 
relative to the non-interacting case, and hence the gravi- 
tational potential is higher, leading to an enhancement of 
growth and a stronger lensing signal. By contrast, when 
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the energy transfer is from dark energy to dark matter, 
< 0, the dark matter density is lower in the past, 
leading to a suppression of growth and a reduced lensing 
signal. In the case r c = 0, the effects of the background 
density evolution are reinforced by the modifications to 
the Hubble friction and source terms in the density per- 
turbation equation (|42p. 

Our conclusions are subject to some caveats: 

• The enhancement or suppression of growth and 
lensing power relative to the non-interacting case 
is tied to the assumption that the current val- 
ues of fl c o and <5 c o are equal in the two cases. 
This gives an accurate picture of how the effects 
of the interaction operate - but we cannot con- 
clude that there is relative enhancement or sup- 
pression, since in reality the current values of Q c q 
and S c o will not typically be equal. The best-fit val- 
ues, based on a range of observations, are likely to 
be different, as shown for other interaction models 
in 0, IIJM [13, H, [U S3 and for the 1^=0 
model in [361 ]. 

• We have followed the standard practice of impos- 
ing the Newtonian approximation on sub-Hubble 
scales: for scales close to the Hubble radius, correc- 
tion terms would start to appear, and numerical in- 
tegration of the exact perturbation equations (|16p- 
(TT8"]) would be necessary. 

• The weak lensing signal that we have computed 
is strictly only valid on linear scales. On nonlin- 



ear scales, the effects of the interaction would need 
to be computed (see, e.g. 0H |H, Q EH) in 
order to construct an accurate lensing signal, and 
this could involve N-body simulations for interact- 
ing dark matter [H, H3, HU . 

Finally, we note that the integrated Sachs- Wolfe effect 
will also carry an imprint of the dark sector interaction, 
and this is investigated for the = model in (3(| H^] , 
and for a different interaction model in 15 Ol . 
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